Naturally inflating on steep potentials through electromagnetic dissipation 
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In models of natural inflation, the inflaton is an axion-like particle. Unfortunately, axion potentials 
in UV-complete theories appear to be too steep to drive inflation. We show that, even for a steep 
potential, natural inflation can occur if the coupling between axion and gauge flelds is taken into 
account. Due to this coupling, quanta of the gauge fleld are produced by the rolling of the axion. If 
the coupling is large enough, such a dissipative effect slows down the axion, leading to inflation even 
for a steep potential. The spectrum of perturbations is quasi-scale invariant, but in the simplest 
construction its amplitude is larger than 10~^. We discuss a possible way out of this problem. 

PACS numbers: 98.80.Cq 



Axion-like particles are the simplest spin-zero degrees 
of freedom with a nontrivial radiatively stable potential. 
Moreover, they are abundant in string theory. As a con- 
sequence, axions provide excellent inflaton candidates in 
a UV complete theory that includes gravity. Inflation 
driven by axions was proposed as early as in 1990 as 
natural inflation The axion potential is radiatively 
stable thanks to a (broken) shift symmetry, and has the 
form = A'*[cos(<l'//) -I- 1], where / is the axion con- 

stant. Neglecting all interactions of $ apart from those 
in and those with gravity, the condition for infla- 

tion is that is flat in units of the Planck scale (i.e., 

|V'| < V/Mp, \V"\ < V/Ml) for a sufficiently wide 
range of In the case of the axion, these conditions 
are equivalent to / >• Mp. Unfortunately, string theory 
appears not to allow such large values of / f^. Moreover, 
/ < Mp appears also as a consequence of the "gravity as 
the weakest force" conjecture of Q. 

In this paper, we show that natural inflation can be 
realized also for a steep potential. Our mechanism relies 
on the coupling of the inflaton to gauge fields through 
the operator ^Ffj^^F'^''. As $ rolls down its potential, 
it provides a time-dependent background for the gauge 
fleld whose vacuum fluctuations are thus amplifled into 
physical excitations. This production of quanta of gauge 
field occurs at the expenses of the kinetic energy of the 
inflaton, slowing it down. If the coupling between $ and 
F^i, is strong enough, such a dissipation effect can allow 
to obtain a sufficiently long period of infiation even if / <C 
Mp. This way, the question of finding natural infiation in 
string theory is formulated in a new way: is there, among 
the many axions of string theory [3] , one whose coupling 
to the gauge fields is large enough? We will argue in 
Appendix B that the answer to this question might be 
positive. 



I. GENERATION OF THE GAUGE FIELD 



We consider natural infiation with a p seudoscalar infla- 
ton $ coupled to a gauge fleld [l^. The Lagrangian 



density of the system is given by 

1 



C 



4f ^^t^'^^ 



(1) 

where V{<^) = A" [1 + cos($//)] with / < Mp 17]. The 
parameter a is a dimensionless measure of the coupling 
of $ to the gauge fleld. 

The equations of motion for the system are 
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where H = a'{T)/a'^{T) and where the prime denotes 
differentiation with respect to the conformal time r. The 
Bianchi identities read B' + 2aHB + y x E — and 
V ■ B = 0. Since the inflaton is homogeneous, V$ = 
0, we can introduce the vector potential A{t, x), with 
a'^B = V X X a^E = -A'. Then, the equations for A 
read 
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In order to study the generation of the electromagnetic 
fleld induced by the rolling pseudoscalar, we promote the 

classical fleld A{t, x) to an operator A (r, x). We decom- 
pose A into annihilation and creation operators 
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where the helicity vectors e± are defined in such a way 
that k-e± = 0, fcx e± — ^i\k\e±. Then, A± must satisfy 
the equations A'^ + {k"^ ^ ak / f)A± = 0. 

Since we are looking for infiationary solutions, we as- 
sume a (r) ~ —\/{Ht), and d^/dt = <i>o =constant. 
Hence, the equation for A± reads 



d^A±iT, k) 
dr"^ 
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A±{T,k) = 0, (5) 
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where we have defined 
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We will be interested in the case ^ > O {1). 

Depending on the sign of ^, one of the two solutions 
or A_ in ([SJ will develop an instability. In the following 
analysis we will assume without loss of generality that 
a > and 6 > (which implies F'($) < 0) so that 

e>o. 

The solution that reduces to positive frequency 
for \k\T -oo is A±{t, k) = [iFo(±^, -fcr) + 

Go(±f, —kT)]/\/2k, where Fq and Gq are the regular 
and irregular Coulomb wave functions. The mode is 
rapidly amplified: when the second term in brackets in 
eq. ([5]) dominates over the first one, |A;t| ^ 2^, A+ is 
approximated by 
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A+ is thus amplified by a factor e^^. On the other hand, 
the modes A_ are not amplified by the rolling infiaton, 
and from now on we will ignore them. 



II. THE SLOW ROLL SOLUTION. 

We can now estimate the backreaction of the gauge 
field on the infiaton, that is described by the term on 
the right hand side of eq. © (note that the backreaction 
of the produced gauge field on the infiaton was studied, 
in a model with different couplings, in Q). Using the 
decomposition of A described above, we get 
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that can be calculated by using the approximate expres- 
sion ([7]) and v4_ ~ 0. Cutting off the integral at kc ~ 
2^i7a(r) and approximating dA^/dr ~ \/2£^ k aH A^, 
we obtain 
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Since we assume ^ >, 1, we will send the upper limit 
of integration in the above equation to infinity, and we 
denote by X = 1\ j (2^'^ i:"^) ~ 2.4 x 10^'' the resulting 
numerical value of the quantity in brackets. 

We then plug (E ■ B) into eq. ^ that, in physical time 
t. now reads 



^2$ xa (H 
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Since we are interested in finding infiationary solutions 
where slow roll is supported by the dissipation into elec- 
tromagnetic modes, we assume that both $ and 3 iJ <I> 



are negligible with respect to V'{^). 
approximate solution of eq. (|10p is 
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where we have assumed 3 Af|> i/^ = ^^'^ j^y {^) + ^{E'^ + 

B^) ~ V{^) (we will check below the regime of validity 
of these assumptions) . Given the logarithmic dependence 
on U($), ^ will never be larger than 0(10). Indeed, un- 
less $ is very close to an extremum of V (and with a not 
exponentially large or small) then £, ^ ^ log [Mp /A] . If 
A ~ 10^ GeV, so that the reheating temperature is much 
smaller than 10® GeV (see below for the estimate of the 
reheating temperature) and overproduction of gravitinos 
is avoided, then ^ ~ 20. 

We can now explore the part of the parameter space 
that leads to infiation. Constraints derive from the fol- 
lowing requirements: 

(i) the Hubble parameter. We first want to approx- 
imate ~ U($)/3Afp, the same relation that holds 
in standard slow roll infiation. This requires that both 
{E^ + B^) and $2 be negligible with respect to U($). 
With the same technique used to estimate (E-B), we 
obtain 
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where we 



have used the slow-roll equation 
/|y'($)|. Eq. (Hg) shows that 
for a 3> ^ the energy in the electromagnetic field can 
be neglected with respect to the energy in the infiaton 
unless we are close to the bottom of the potential. 
Indeed, by approximating ^('I') oc $2 near its minimum, 
we see that when $ < $rh = C//q^ the energy in elec- 
tromagnetic modes cannot be neglected any more. This 
is the point where reheating begins. The energy density 
at reheating is ~ A*$^jj//2, so that the reheating 
temperature will be of the order of A yj^ja. 

Next, let us analyze the condition $2/2 ^ V . Using 
$ = 2fH^/a, we obtain 
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that shows that the kinetic energy of the infiaton can 
be neglected with respect to the potential energy for 
(C/a) if/Mp) < 1. Since f <Mphy assumption, a > f 
will be sufficient to satisfy this condition too; 

(ii) acceleration. In order to make sure that our so- 
lution actually corresponds to an inflating Universe, we 
compute the slow roll parameter e = —H/H^. Using the 
equations of motion ([2]), we obtain 
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where isotropy of the background imphes (V • x B)) = 
0. By inserting into eq. the values of 4>^, {E'^ + 
and found above, we derive the expression 



field as $ = $0(7") + (/'(r, f), then C = H (j)/^o- In order 
to compute the power spectrum of ^, we must therefore 
compute the correlators of 0. (f) obeys the equation 
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By comparing the above equation with eq. (fT2|) and (|T3|) 
we see that, as long as the conditions va (i) are satisfied, 
we will have e < 1 and the Universe will be inflating; 

(Hi) neglecting terms in $ and 3 H ^ in eq. I110\}. The 
following conditions must be satisfied: 



(a) 



V 



e fv/v 

2 a Ml 



«1 , 
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Since /F/^' = 0{p) and V^V^''/^'^ = 0(1) unless we 
are close to an extremum of the potential, then £^^\ 
guarantees that both (a) and (b) hold; 

(iv) number of efoldings. The strongest constraint 
comes by requiring that inflation lasts for long enough. 
The number of efoldings is given by 
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a^S[E-B], (18) 



where the fluctuation d[E ■ B]{t,x) receives two contri- 
butions. Besides the intrinsic inhomogeneities m. E ■ B 
(that would be present also for (j) — Q), a second compo- 
nent comes from the fact that {E ■ B) depends on $. As 
a consequence, if $ is replaced by $-!-(/', then {E ■ B) will 
go to {E ■ B) d{E ■ B) We therefore write 

5[E-B]^[E-B-{E- B)]^^o + ■ <P (19) 

We denote the term in square brackets by Sg g{r, x). 
The second term is estimated by observing that {E ■ B) 
depends on i> through ^ and that d{E-B)/d^ ~ 2n{E-B). 
Using the background equation a{E ■ B) / f ~ V , the 
second term of the right hand side of eq. (fT9)) can be 
written as tt aV (j)/{fH). The Fourier transform of the 
perturbation will then obey the (operator) equation 
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Since the range of variation of $ is bounded by |$/ — 
\ the above equation implies that a ^ 2^iVe/7r. 

Hence, for ^ ~ 20 we need a ~ 600 to obtain 45 efoldings 
of inflation [H]. 

To sum up, natural inflation with electromagnetic dis- 
sipation will last for efoldings ii a ^ 2 ^ Ne/n, where 
£_ - {2/tt) log(10a^i/'*Mp/A). There are no limits on 
the scale of inflation A (apart from the obvious require- 
ment A >, GeV to allow for Big Bang Nucleosynthesis 
after reheating). Reheating will occur when $ = 'I'rh, 
with 
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III. PERTURBATIONS 



Perturbations are usually generated as the quantum 
fluctuations of the inflaton are amplified by the evolving 
background. This is different from our scenario, where 
inhomogeneities in $ are sourced classically by those in 
the electromagnetic field. The situation is analogous to 
the one studied in [1, 0], and we will use similar tech- 
niques to analyze it. 

The curvature perturbation ^ on a uniform energy den- 
sity final hypersurface is related to the perturbation of 
the number of efoldings by C = SN = N{x) — N, where 
N is the number of efoldings in the homogeneous back- 
ground. If we write the perturbed value of the inflaton 
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Denoting by G (r, r') the retarded Green function as- 
sociated to the differential operator acting on in the 
equation above, the two-point function of the inflaton in 
momentum space reads 
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where we use the notation a' = a (r'), a" = a (t"). 

Therefore, we must compute the two-point correlator 
of g, find the Green function associated to the homo- 
geneous part of eq. (j20l) . and compute the integrals in 
eq. (PT|) . The details of this derivation are presented in 
Appendix A. The final result is 
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where (see appendix A) 



(22) 



v+ ~ T^aV'lif H^) cx a Ml/f > 1 , 

v_ ~ V"f/{T:aV') cx 1/a < 1 . (23) 

The curvature perturbation is = 
iJ^ ((/)(/))/ [2 TT^ $0 J(p + p')]. For generality, we 
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extend the result to the case where the theory contains 
Af gauge fields. It is straighforward to see that for 
Af 1 the constraints on the parameters found above do 
not change. However, the different contributions to the 
two point function of add incoherently, so that V(; 
is suppressed by a factor l/f\f. Taking this suppression 
into account, and using a(_ff/^)^ e^'^^ = f\V'\/I, we 
finally obtain 
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The spectral index of the scalar perturbations is 
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While the sign of V does not change during infiation, 
V" crosses zero. As a consequence, the spectrum can be 
either red or blue depending on the value of $ at the time 
the relevant scale exited the horizon. 

The amplitude V^^ - 0.05/A/'C^ matches the COBE 
normalization 2.5 x 10~^ only for large values of A/", since 
^ = 0(10). In particular, if we assume ^ ~ 20, then we 
need A/" ~ 5 x 10^ to obtain perturbations with the ob- 
served amplitude. Such a large number of gauge fields 
can be obtained for instance if the theory contains M 
branes, each with its own U{1) gauge field. Alterna- 
tively, one can obtain M gauge fields by considering a 
gauge group SU{VJ7) (for instance, ref. [l3] considers 
groups as large as S'J7(520)), that might be obtained on 
a stack of ^/Af ~ 200 branes. This option has the advan- 
tage that the different gauge fields have automatically the 
same coupling a to the infiaton. Note that in the case of 
nonabelian groups, one should in principle take into ac- 
count the self interaction of the gauge fields. However, as 
long as the gauge self-coupling is weak, its effects appear 
only at higher order, and can be consistently neglected. 



IV. CONCLUSIONS AND FUTURE 
DIRECTIONS 

We have shown that it is possible to realize infiation 
on a steep axionic potential, provided the infiaton has 
a sufficiently strong coupling (of the order of ^ (10^ — 
10'^)//) to a gauge field. Remarkably, this scenario can 
accommodate inflation at any energy scale and for any 
value of / Mp. Unfortunately, the simplest version 
of the scenario gives an exceedingly large amplitude of 
scalar perturbations. With a sufficiently large number 
of gauge fields it is possible to reduce such amplitude 



J 



to the observed value. It would be interesting to see 
whether such a suppression could be achieved by other 
mechanisms. 

We still need to show that a as large as 10^ — 10'^ 
can be realized. In Appendix B we show two examples 
(one directly derived from a string construction, the other 
based on extra dimensions) where such large values of the 
coupling can be obtained. 

Our mechanism is in spirit similar to that at work in 
warm fT2l| and trapped inflation [l^ , that also use dis- 
sipation to realize slow roll. The main difference from 
these models is that our infiaton is coupled to a deriva- 
tive to the produced field. This allows to achieve a sta- 
tionary dissipative process without relying on the more 
complicated structures invoked in those models. 

A number of details of this scenario still need to be 
explored. First, since perturbations are sourced by the 
inhomogeneities in the gauge field, it is important to 
study nongaussianities in this models. Future work 
also involves the generation of gravitational waves, and 
it will be interesting to study whether a consistency 
relation similar to that of standard slow roll inflation 
holds also in this case. One flnal question concerns 
parity violation. Since our electromagnetic field is 
maximally parity violating [l^ . the gravitational waves 
produced by the electromagnetic modes could generate 
a nonvanishing (TB) correlation [l^ in the Cosmic 
Microwave Background. 

Acknowledgments We thank Nemanja Kaloper, David 
Langlois, Anthony Lewis and Alessandro Tomasiello for 
discussions. This work is partially supported by the NSF 
grant PHY-0555304. 



Appendix A: Calculation of the power spectrum 

As discussed in section HI, the curvature perturbation 
power spectrum is related the two point function of the 
perturbations of the scalar 0, eq. (j2ip . In this appendix 
we derive the power spectrum starting from that ex- 
pression. In order to perform this calculation, we will 
first need to compute the two point correlator of Sg 
and the propagator G{t, t'). 



1. The two point function of 5g g 

By using the definitions and the properties of section 
3 and assuming A_ ^r, A;^ ~ 0, the correlator is given by 
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In principle, the above expression should be renormal- 
ized. However, the only part that matters to the gener- 
ation of inhomogeneities in (f> corresponds to the wave- 
lengths larger than (2^ H)~^, where the electromagnetic 
field has large occupation numbers and can be treated as 
a classical source. In this regime, the function ^+(t, k) 
is well approximated by eq. ([?])• As a consequence we 
can write the above equation as 



as / ;S Mp/^, is always satisfied if p <C 2^aH, that is 
the necessary condition for the validity of eq. ([7]). 

Once we have taken into account these two approx- 
imations, we can find the appropriate Green function, 
obtained by solving 
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where we have defined the function a (r', t") via 2/\/a — 
l/v^ -I- 1/vV, with a' = air'), a" = air"). Note that 
the integral in the above equation (|A2p should be cut- 
off at \p + k\ < 2^HMin{a', a"}. However, using a 
reasoning analogous to the one made after eq. ([§]), we will 
send the integration limit in k to infinity. 

After a change of integration variable, we finally write 
the correlator as 
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G'(t, r') = 5(r-T') . (A5) 
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where the function C(k) (after directing p along the z 
direction) reads 
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where z is the unit vector in the z direction. 



2. The Green function 

In principle, the Green function for equation (|20|) can 
be computed exactly (after assuming $o —constant). We 
can however get simpler analytical expressions by limit- 
ing ourselves to the case of the cosine potential T^($) oc 
1-f cos($//). In this case, V ($o) ~ l^($o)//, V" (*o) - 
F($o)//^ -ff^ ^ V{^q)IMI and a > 1 while / < Mp. 
This allows to see that the coefficient of d(t>/dT in eq. ((20)) 
is much larger than one. Moreover, we can neglect the 



term p^ i n th e coefficient of cj). This is possible for 
\p\ <^ a \/\V"\ ~ {Mp/ f) aH , a condition that, as long 



with G(r', t') = 0, (9G/9r)(T', r') = 1. The solution is 
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The second term under the square root in eq. (|A7|) scales 
as {f /aMpY and is much smaller than one. As a con- 
sequence we have v+ ~ naV / {f H'^) oc aMp/p > 1 
whereas ~ ^"//(TraT^') cx l/a < 1. 



3. The spectrum 

We are interested in the spectrum at p <C aH . As a 
consequence, we can neglect the term {t/t'Y+ in the ex- 
pression of the Green function, that goes rapidly to zero. 
We then collect the results of the above subsections, we 
use a = — 1 / Ht and define the new integration variables 

w' = ~ (2^^ \p\ t')"\ w" = - (2^^ IpI r")~\ There- 
fore, the two point function reads 
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where we have defined 2/\^ = l/\/u7+ 1/Vw". We see 
now that, as long as <C 1 the spectrum of perturba- 
tions in the infiaton is quasi-scale invariant. To find the 
normalization, we send p — > in the limits of integration. 

The integral can then be computed by defining the 
new integration variables x' — x" — w"^^/'* and 

going to "polar coordinates" x' = pcosO, x" — ps'md. 
We set = to simplify the resulting expressions. We 
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use the expression (jA4l) for the function C. The integrals 
in 9 and p can be now computed explicitly 
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so that the fields 3> and ^ are canonically normalized and 
the interactions between the axions and the F^n and G^,y 
read 
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where the integral in d'^q can be evaluated numerically 
to ~ 3.5 X lO-'^. 
The final result is 

2i/_ 
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where the numerical factor 7~2.1xl0 ^ 

We can finally find the curvature perturbation Vq = 
H^(b<j)) /[2Tr^^ 6^^^ {p + p% Using a(i?/^)4 e2-« = 

fV /I, we obtain 
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(All) 

If the theory contains Af gauge fields, then the different 
contributions to the two point function of <5^.^ will sum 
incoherently, leading to a suppression by a factor of J\f of 
•Pc, i.e.. 
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Appendix B: Obtaining a = 0(100) 

We present here two examples where one can obtain a 
value of a large enough to allow our mechanism to work. 
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The field <& can be ignored, and if the G^i, sec- 
tor gets strongly coupled at energy A, then a po- 
tential ~ A'' cos(<i>//off) is generated. By identifying 
(/I + /i ) / (/I ^ /i ) = a we can obtain a 1 by fine- 
tuning /i ~ /a (1 - 1/a). 



2. Extra dimensions 

A second possibility involves extra dimensions. To 
fix ideas, let us consider a flat, five dimensional bulk 
compactified with radius R. A single Peccei-Quinn field 
jy — ^3/2 gi e jjygg |;jjg bulk, where we have already 
fixed the absolute value of rj to v^/"^ and we will neglect 
Ty's radial mode. Assuming that the kinetic term of rj 
is canonically normalized, we have that the lagrangian 
for a reads — / d^xdyv^ {80}^. Let us now consider two 
gauge fields Aj^ and A^, both living in the bulk. Both 
gauge fields will have an axionic coupling to 9, but while 
the coupling to Aj^ will be in the bulk, the coupling to 
Aj^ will be localized on a brane at y = 0. Therefore, the 
relevant part of the lagrangian reads 
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1. Two axions, two gauge fields 

The first example we consider is based on a lagrangian 
presented already in [ll| , that studied axion phenomenol- 
ogy in the dimensionally reduced x superstring. 
The low energy theory contains a model independent ax- 
ion $1 and a model dependent one, $2 
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where i^^^ and C^^ denote the field strength of two dif- 
ferent gauge groups. 
We then redefine 
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/2 



/i_ 

fi 



/2 



fi 



(B2) 



where M is some mass scale that we assume to be of the 
order of the cutoff scale of the theory. We obtain a four 
dimensional effective action by integrating on the extra 
dimensions, and canonically normalizing both the axion 
and the gauge fields. The resulting lagrangian for the 
zero modes reads 



(dOf 



-F' F 
4 ^ 



1 flU 



1 



(B5) 



p2 p2^.v 



Let us now assume that the second gauge group con- 
denses at some scale A. Then a potential ter m ^ 
A^ cos{9/f) will be generated, with / ~ M R y/v^ R. 
The coupling of 9 to the first gauge group will thus have 
the form M R{e/ f) F^^ F^f"" so that a ^ MR can be 
of the order of 10^ provided the extra dimension is large 
enough. 
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